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Abstract 

We construct an SO (5) symmetric electron model on a two-chain ladder with 
purely local interactions on a rung. The ground state phase diagram of this 
model is determined in the strong-coupling limit. The relationship between 
the spin-gap magnon mode of the spin-gap insulator and the tt resonance 
mode of the d-wave pairing phase is discussed. We also present the exact 
ground state for an 50(5) superspin model. 



I. INTRODUCTION 



The cuperate two-leg ladder materials are characterized by strong electronic correlation 
and a variety of competing ground states |l| . Numerical and analytic calculations for Hub- 
bard and t — J models || have shown that at half-filling, these two-leg ladder models 
exhibit a spin-gap insulating phase and that when the system is initially doped, rf x 2_ 3/ 2-like 
pairing and CDW correlation can become dominant. At higher doping the system is ex- 
pected to behave as a one-dimensional Luttinger liquid M. Given the strongly interactive 
nature of these systems and the delicate balance between their competing ground states, 
one would like to have a more general framework for determining their properties. Recently 
it was suggested that one could capture some of the basic low energy physics with mod- 
els having an 50(5) symmetry |5]||. As we will discuss, there is in fact a natural way to 
construct an 50(5) symmetric model for a two- leg ladder which has only local interactions 
on a rung. In this work, we will show how to construct such a Hamiltonian and discuss its 
strong-coupling phase diagram and low lying collective excitations. 

This work is partially motivated to use the ladder system as a theoretical laboratory to 
check some ideas of the 5*0(5) theory. Recently, two dimensional lattice electronic models 
with exact 5*0(5) symmetry have been constructed by three groups independently 0-0. 
However, these models all involve long range interactions. The 50(5) symmetric ladder 
models constructed in this work involve only local interactions on the rung and are therefore 
much easier to visualize. There has been considerable progress in numerically checking the 
approximate 50(5) symmetry in the 2D t — J and Hubbard model ||TO| , pT| . The locally 50(5) 
symmetric ladder models constructed in this work are simple to implement numerically. By 
systematically varying the parameter away from the 50(5) symmetric point, one can trace 
the evolution of the 50(5) multiplet structure |Tl| and get a better sense of the nature 
of the approximate 50(5) symmetry. In this work, we find a continuous quantum phase 
transition from the spin gap Mott insulator phase to the d wave superconducting phase, 
and we show that the spin gap magnon mode of the Mott insulator evolves continuously 
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into the 7r resonance mode [|6||T0l-|T2|l of the superconducting phase. These results also shed 
some light on the nature of the ir resonance mode in the 2D case. Perhaps one of the 
most important questions in the 5*0(5) theory concerns the origin of such a symmetry in 
generic models. Recently, Shelton and Senechal ]13[ studied the problem of two coupled ID 
Tomonaga-Luttinger chains and concluded that approximate 5*0(5) symmetry can emerge 
in the low energy limit of this model. Balents, Fisher and Lin jl4| have used a weak 
coupling RG method combined with abelian bosonization to show that a generic ladder 
model at half-filling flows to a manifold with 50(5) symmetry (and in fact, to various phases 
characterized by higher symmetries). Arrigoni and one of us (WH) additionally included a 
next-nearest-neighbor hopping t', which explicitly breaks the symmetry between the bonding 
and antibonding bands and thus the 50(5) symmetry even in the non- interacting limit ||15|| . 
In this case, the model still flows to an 50(5) symmetric (at least up to order 0(t'/t) 2 ) 
effective action, provided the 50(5) symmetry is redefined by accounting for different single- 
particle renormalization factors for the fermions on the two bands of the ladder system. 
These results indicate that the 50(5) symmetry can emerge as a result of RG flow, and it is 
therefore of interest to study the low-energy physics of fixed point Hamiltonians which have 
exact 50(5) symmetry. 

This paper is structured as follows. The formal construction of an 50(5) symmetric two- 
leg ladder Hamiltonian will be described in Section II. Following this, Section III contains a 
discussion of the ground state phases in the strong coupling limit. The collective modes are 
discussed in Section IV, and Section V contains a summary of our results. 

Before going into the technical details of the subsequent sections, we would like to con- 
clude the introduction by explaining the basic idea. On a given rung of the ladder, there are 
two sites with 16 electronic states, as depicted in Fig. 1. These 16 states can be classified 
into 4 different groups, (the E , Ei, E 2 and E 3 groups in Fig. 1), each transforming as 
irreducible multiplets under 50(5) [[LI]]. If a local Hamiltonian on a given rung is 50(5) 
symmetric, states within a given multiplet must have the same energy. However, simple 
visual inspection shows that the states in the E 2 and E% groups are already degenerate 



for generic interactions respecting spin invariance, particle-hole symmetry with respect to 
half-filling and symmetry under interchanging two sites. Therefore, only one condition is 
required to make the spin triplet magnon state and the two "pair states" in the E\ group 
degenerate. This condition turns out to be J = 4(U + V), relating the onsite interaction U, 
a near-neighbor interaction V, and an spin exchange interaction J on a rung. 

II. CONSTRUCTION OF 50(5) SYMMETRIC LADDER MODELS 

Rabello et al. showed that the crucial step in constructing general exact SO (5) sym- 
metric models is to identify a 4 component fermion operator which transforms according to 
the fundamental spinor irreps of 5*0(5), or its equivalent Sp(A). The important point here is 
that the geometry of the 2-leg ladder makes it natural to group fermion operators on a rung 
into a 4-component spinor, which leads to a Hamiltonian with purely local interactions. In 
the following we will use x, y, .. to denote the position of a rung on a ladder, and c a (x) and 
d a (x) to denote the spin up (a = 1) and spin down (a = —1) fermion destruction operators 
on the upper and lower chain, respectively. Our 50(5) spinor operator is defined by 

* a (x) 

for the even rungs with (— l) x = 1, and 

* a (x) 

for the odd rung with (— l) x = — 1. These spinor operators satisfy the canonical anticom- 
mutation relations: 

{*t(x),* fl (y)} = 6{x-y)6 a p (2.3) 

and 

{*«(*), **(y)} = {*!,(*), *Uv)} = ( 2 - 4 ) 
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( c„(x) ^ 



' d a (x) ' 



(2.1) 



(2.2) 



Using these spinor operators and the Dirac T matrices (for details see appendix A), we can 
construct a 5 dimensional 5*0(5) superspin vector 

n a (x) = U>i{x)Tl^> p {x) (2.5) 

a 10 dimensional SO (5) symmetry generator 

L ab (x) = -U>i{x)T%^ p {x) (2.6) 

and a 1 dimensional 5*0(5) scalar 

p{x) = (2.7) 

on a given rung x. The local commutation relation between these operators are given by 

[L a b, L cd ] = —i(5 ac L bd + 5 bd L ac — 5 ad L bc — 5 cc L ad ) (2.8) 



[Lab, n c ] = -i(5 ac n b - 5 bc n a ) (2.9) 

and 

[L ah ,p} = (2.10) 

The superspin vector n a is related to the AF and SC operators by 

( At + A ) 1 + + 

ni = A '- = -{-iJo-yS + h.c.) (2.11) 

n 2 , 3 ,4 = N x%VtZ = X - {Sac - Sad) (2.12) 
(At -A) i 

n 5 = V ^ 1 = ~-{c ] a y S + h.c.) (2.13) 

where we supressed the spinor index on c a and d a . The symmetry generators L ab are 
expressed in terms of the rung spin S a = \{Sa a c + d^a a d), charge Q = \{Sc + d)d — 2) and 
the 7r Q operators 

rrl = -\ja a a y S (2.14) 
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with 



7rl + 7T x 





-s z 



7U + 7T 2 







(2.15) 



/ 



^ Q l -(nt-n x ) i(7rJ- % ) \(4-* M ) 

The 5*0(5) singlet operator p = \{c*c — cftd + 2) has the physical interpretation of the 
charge- density- wave operator. 

Having exhibited the local 50(5) operator algebra, we are now in a position to construct 
an 50(5) symmetric model. Let us first consider the problem of two sites on a given rung. 
As discussed in the introduction, there are 16 states which can be classified under 50(5) as 
1) a spin singlet state on the rung 



I") = v ^ J |Q) 



(2.16) 



which is also an 50(5) singlet (L a b|f2) = 0), 2) an 50(5) vector quintet n a \Q) which contains 
a triplet of "magnon states" and a doublet consisting of a hole pair and its conjugate, 3) 
two 50(5) spinor quartets \J/ Q |^) and \frf|fi), and finally 4) two additional 50(5) singlets of 
the form ^ a R a ^ ^Q,) and ^i? a/ j^t|fi). The R matrix is a invariant tensor of the 50(5) 
algebra and it is defined in the Appendix A. These states are depicted in Fig. 1. 

Let us first neglect the hopping within the rung, and consider the general spin and charge 
interaction Hamiltonian for the two sites: 



runq 



1 1 



= U (n aT 
+ V(n a -l)(n b -l) + JS a S b 



(2.17) 



In order for such a Hamiltonian to be 50(5) symmetric, we would require that the states 
within each of the four multiplets mentioned above be degenerate. As noted earlier, de- 
generacy within the ^ a \fl), \I/JJQ), multiplets is automatically ensured by spin rotation 
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invariance, the particle- hole symmetry and the symmetry under rung parity (c CT — > do-). The 
^ a Ra/3^ an d ^a-Ra/?^]?!^) states are already 5*0(5) singlets. This leaves only the n a \Q) 
quintet manifold. For the rung Hamiltonian, Eq ( |2 . 1 7| ) , the triplet "magnon states" are seen 
to have energy J/4 — U/2, while the doublet pair states have energy U/2 + V. Therefore, 
these states will be degenerate if 

J = A(U + V) (2.18) 



Under this condition, the Hamiltonian ( |2.17 ) can be cast into the manifestly 5*0(5) sym- 
metric form: 

Hrun 9 = jE^+(| + - 2 ) 2 ( 2 - 19 ) 

a<b " 

up to an additive constant 3 J/4 + U/2. For this 50(5) symmetric form of the interaction, 
the above mentioned four multiplets have the energies E = —\U — 3V, E\ = U/2 + V, 
E 2 = and E 3 — U/2 — V respectively as indicated in Fig. 1. 

In the manifestly 50(5) symmetric Hamiltonian fl2.19|) , the interactions are expressed in 



terms of tensor and scalar interactions. The readers may wonder why the vector interactions 
are missing. The reason is that there is a Fierz identity discussed in Appendix A, which 
relates the three channels of interactions, and only two of them are mutually independent. 

Let us now consider the effect of the hopping within the rung. It is easy to see that the 
hopping term can be expressed in the manifestly 50(5) symmetric form: 

H tx = -2t ± (cld a + h.c) 

= t ± (ty a R a ^p + h.c.) (2.20) 

This hopping term can split the degeneracy within the E 2 and E 3 manifold. The E 2 manifold 
splits into anti-bonding and bonding states (\l/ a ± (R^) a )\Q), with energies E% = 2t± and 
E\ = —2t± respectively. H t± can in principle cause mixing between the Eq and the two 
states in the E 3 manifolds. Since all three states are 50(5) singlets, their mixing does not 
violate 50(5) symmetry. Within degenerate state perturbation theory, the two E 3 states 



also form anti-bonding and bonding combinations (^ a R a p^p i^^)^), with energies 
E^ — U/2 — V — t\/\2U + V\ and E\ — U/2 — V + t\/\2U + V\ respectively. 

So far the 50(5) symmetry is only realized on the two sites of a rung. The more non- 
trivial question is how the symmetry is realized when the hopping tn in the ladder direction 
is included. Remarkably, this hopping can also be expressed in manifestly 50(5) invariant 
form, 

H H = — £ (ct(x)c a (y) + 4(x)d a (y) + h.c) 

<x,y> 

= 2t|, (*a(x)R a ^f3(y) + h.c.) (2.21) 

<x,y> 

It is important to point out that the alternating definition of the 5*0(5) spinors on the even 
( p.l|) and odd (|2.2|) rungs makes it possible to express H t in manifestly 50(5) invariant form. 



This alternating definition was suggested to us by S. Rabello. Without such an alternating 
definition, the resulting Hamiltonian is still 50(5) symmetric, but the symmetry is not 
manifest. Our final 50(5) symmetric ladder Hamiltonian H is given by the sum of ( [2.19P , 
( p.20| ) and ( |2.21| ). For this Hamiltonian, all ten 50(5) generators 



L ab = J2 L *b{x) (2-22) 

X 

are exactly conserved, 

[H,L ab ] = (2.23) 



Notice that because of our alternating definitions ( |2.1|) and ( |2.2| ) of the fermion operators, 
the 7r a operators have momentum 7r along the ladder, while the total charge and total spin 
operators are uniform. In the presence of a chemical potential term, = — 2/iL 15 , the n a 
operators are exact eigen-operators 

'h + H^ttI] =-2 M 7rt (2.24) 

so that the total Casimir charge of the ladder C = J2 x ,a<b L 2 a b ls conserved. 

[H + H II ,C} = (2.25) 



Therefore, all states of the doped ladder are still labeled by their SO (5) quantum numbers. 

The 5*0(5) symmetric ladder model we presented so far has only local interactions on the 
rungs. Obviously, one can generalize the model by including interaction between different 
rungs, for example one could write down SO (5) invariant interactions having the form 

V x {x - y)n a (x)n a (y) + J2 v i( x - y)L ab (x)L ab (y) 

x,y x,y 

+ £ V (x - y)(p(x) - 2){p{y) - 2) (2.26) 

x,y 

Here we shall restrict ourselves only to the analysis of models with local rung interactions, 
and defer the general analysis to future works. It is plausible that in the strong coupling 
limit, the local rung interaction dominates the physics. 

III. STRONG COUPLING PHASE DIAGRAM 

In this section we discuss the phase diagram of the SO (5) ladder Hamiltonian in the 
strong coupling limit. Setting J = A(U + V), the energies of the different rung manifolds 
are listed in Fig. 1. One can divide up the U-V plane according to regions in which a given 
manifold lies lowest in energy and Fig. 2 shows such a plot. In the strong coupling regime, 
one can study a given sector of the U-V plane, using the virtual hopping processes due to 
H t± and H t ,, to resolve the degeneracies and determine the dynamics of the low lying excited 
states. In the following we examine the three different regions E , E 3 , and E\ shown in the 
U-V phase diagram of Fig. 2. 

A. The E Spin-Gap d-Wave Phase 

In the region E bounded by V = —2U for negative values of U and V = —U for positive 
values of U, the singlet rung state \Q) (see Fig. 3a) is lowest in energy and the system is 
expected to be in a spin-gap insulating ground state. In the strong-coupling limit this state 
is simply a product of rung singlets. The spin-gap corresponds to the energy to create a 
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magnon triplet, as illustrated in Fig. 3b, in which a rung singlet is replaced by a magnon 
triplet from the E\ manifold. This costs an energy 



^s 9 = ( ¥ + v 



7U \ 

— + 3V) = A(U + V) = J (3.1) 



If we were to add two holes to this phase, the lowest energy state occurs when the two 
holes are placed on the same rung as illustrated in Fig. 3c. In this case, the excitation 
energy is again A sg = J, as expected for an SO (5) symmetric system. Alternatively, one 
could imagine adding two holes by placing each one on a separate rung creating two E 2 = 
states. However, this would cost an energy 2(| U+3V), because of the two singlet rung states 
that are destroyed, and this is a larger cost in energy than placing the holes on the same 
rung throughout the entire E region of Fig. 2. Thus the doped holes will form rung pairs 
and we expect that the doped system will exhibit power law pairing and CDW correlations 
in the E region. If one defines bonding and antibonding rung orbitals 

bUx) ^M±m 4{t)= m^m (3 . 2) 

then the singlet rung state has the form 

c\d\ — cldf\ . . ( btb\ — a\a\ , , , 

i I | ) = T i T 1 |0) (3.3) 



y/2 ] \ V2 

Thus in the E'o-region, the hole pairs go into a "c?-wave" like state [Q] in which the amplitudes 
of the singlet pair in the bonding and antibonding orbitals have opposite signs. 

Both the magnon and the hole pair can propagate coherently along the ladder leading 
to an energy dispersion in q x . In strong coupling we can calculate their dispersion relations 
to second order in tu as follows. A magnon excitation on rung x can hop to rung x + 1 by 
going through an E\ intermediate state. If 1^) is a state with the magnon on site x, then 
the second order virtual hopping process has a matrix element 

and the magnon dispersion is 
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2 



3U + 2V 



cosq x 



(3.5) 



There is also a £jj shift in the zero point energy when a magnon is created. Taking these 
virtual processes into account gives the complete magnon dispersion to second order in tfr. 



It is straight forward to carry out a similar calculation for the hole pair dispersion and 
one finds that only the sign of the cosq^ term in Eq. (|3.6|) is changed. Thus, as expected 
for an SO(5) symmetric ladder, the magnon dispersion about q x = it is identical to the hole 
pair dispersion about q x = 0. 

If two hole pairs (or two magnons) are added, one finds that they repel each other 
if they are located on neighboring rungs. This repulsion simply reflects the reduction in 
the zero point fluctuation energy from the background fermions when the two hole pairs 
are adjacent. This gives rise to an effective near neighbor repulsion V* = 4t|/(3C7 + 2V). 
Therefore, in the Eq region of the U — V phase diagram the doped system behaves as a dilute 
one-dimensional boson gas with a repulsive near-neighbor interaction, and is expected show 
quasi-long-ranged-order in the superconducting correlation function. The transition from 
the spin gap Mott insulator to the superconductor occurs when the chemical potential is 
given by 



In the i? 3 -region of the U-V strong coupling phase diagram (Fig. 2) the singlet states 
^ttRai3^p\Q) and a R a pty are coupled by second order kinetic energy processes. In this 
regime, the problem can be mapped onto an effective Ising model in a magnetic field. We 




(3.6) 



2\fi\ = 2fi c = A 



(3.7) 



and such a transition is expected to be second order. 



B. The £3 Spin-gap s-wave and CDW Phase 
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identify ^ a R a /3^\Cl) with the Ising state \a z (x) = 1) and \& a Ra0^ with the Ising state 
\a z (x) = — 1). On a rung we have 

<**<*> = ^e±b;W) = -D = lihv < 3 - 8 ' 

and between two near neighbor rungs 

In the i?3-region, the energy denominators in the expressions are negative so that the t\ 
term favors the formation of the s-wave like rung singlet 

+ + , , , , b\b\ + ala\ , , 

¥ a R a ^l\n) + t> a R a ^f3\n) = T 1°) ( 3 - 10 ) 

Here and aj. are the bonding and antibonding creation operators of Eq (|3.2|) . On the other 
hand, the t| process favors a staggered charge density wave state. Combining equations ( p.7|) 
and ( |3.8| ) we can write an effective Ising-like Hamiltonian for the ^-region in the form 

^3 = E (-ho*(x) + K 3 (a z x+1 a z {x) - 1)) (3.11) 

X 

with h = t\/\2U + V\ and K 3 = 2t\/\U/2 - V\. The ground state of H 3 is known to have 
an Ising-like phase transition for h = K 3 . For h < K 3 , the half-filled 50(5) ladder will be 
in a CDW phase corresponding to one of the two degenerate states illustrated in Fig. 4a 
and Fig. 4b. For h > K 3 , the system will be disordered. In this region, the half-filled 
5*0(5) ladder will be in a spin-gap insulating phase. For tu = t±, the h = K 3 dividing line 
corresponds to V — — jU. When holes are doped into the disordered region, they will tend 
to go onto a rung forming s-wave like, Eq ( |3.10|) , pairing correlations. 

C. The Ei 50(5) Superspin Phase 

Here the Hilbert space per rung is restricted to the "superspin" quintet manifold 
n a {x)\Q). In this case, each rung is either occupied by a triplet magnon or a doublet "pair" 
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state. The effective Hamiltonian in the quintet manifold is easily determined using second 
order perturbation theory: 

Hi = Kt ]T L ab {x)L ab {y) (3.12) 

<x,y> 

where Ki = t?J\U/2 + V\. This model can be viewed as the SO(5) generalization of the 
spin one Heisenberg chain. Therefore, we would expect to find many properties simply from 
this analogy, for example, a ground state with a finite excitation gap, and short ranged 
correlations, etc. One useful model of the spin one Heisenberg chain is the AKLT ]16|] model 



for which an exact ground state is known. In this section we shall construct an SO (5) 
generalization of the AKLT model and present its exact ground state. 

We begin by considering two neighboring rungs x and y. The wave function for the two 
superspins defined on the two rungs can be decomposed as 

5x5 = 1 + 10 + 14 (3.13) 

i.e. the product wave function can transform like an SO(5) singlet, an SO(5) antisymmetric 
tensor or an SO(5) symmetric traceless tensor. Therefore, we can defined a complete set of 
bond projection operators P\(xy), P\o(xy) and Pi4,(xy) onto these subspaces, satisfying: 

1 = P 1 (xy) + P w (xy) + P 14 (xy) (3.14) 
The 5*0(5) generalization of the AKLT model is then given by 

H X = 2K X p i^y) ( 3 - 15 ) 

<x,y> 

Let us first see how this Hamiltonian can be expressed in terms of the SO (5) superspin 
exchange operators. We start by defining 

£ab = L ab (x) + L ab (y) (3.16) 

where C ab measures the total 50(5) generator on the bond < xy >. Squaring this equation 
and noticing that the Casimir charge O5 = J2 a <b L 2 a b( x ) f° r the 5 irreps on a given rung is 
O5 = 4, we obtain 
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EM^{y) = 5E^- 4 ( 3 - 17 ) 

a<b a<b 

The operator J2a<b £-ab * s the total Casimir charge for the bond < xy >, therefore, it can be 
expressed as 

£lb = CiPi(xy) + C 10 P 10 (xy) + C 14 P 14 (xy) (3.18) 

a<b 

where C\ = 0, Cw = 6 and C\ 4 = 10 are the Casimir charge for the 1, 10 and 14 irreps 
respectively. Therefore, we obtain 

£ L ab {x)L ab (y) = -AP^xy) - P 10 (xy) + P 14 (xy) (3.19) 

Squaring this equation again gives 

(]T L ab {x)L ab {y)) 2 = 16P 1 (xy) + P 10 (xy) + P 14 (xy) (3.20) 

where we used the property of the projection operators P%Pj = bijPi- Equations ( p,14[ ), 
( |3.19j ) and ( |3.20| ) finally allows us to express the P\ 4 {xy) operator as 



Pi4*v) = ^(L ab (x)L ab (y)) 2 + ^(L ab (x)L ab (y)) + ~ (3.21) 



Inserting this equation into ( |3.15|) then gives the desired expression for the SO (5) general- 



ization of the AKLT Hamiltonian. 

The Hamiltonian ( |3.15| ) has an exact ground state. In general, the ground state can be 



expressed as 

1^0)= i>a( a u-; a N)n ai {xi)--n aN (x N )\ty (3.22) 



ai,..,ajv 



for a ladder with N rungs. ^o(ai, .., ajy) is the corresponding ground state wave function in 
the superspin vector basis. It is easy to see that the exact ground state wave function for 
the SO (5) AKLT Hamiltonian ( ETT5 ) for a period ladder is given by 



Mai,a 2 ,..,a N ) =Tr{r^T a2 ..T aN ) (3.23) 
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This follows from the following property of the Dirac T matrices (for more details, see 
appendix A): 

T a T b = 26 ab + 2iT ab (3 34) 

From this equation we see that the product of two T matrices involves no symmetric traceless 
components. Therefore, the wave function ^o(«i ; ••, a-i, 04+1, •-, o>n) viewed as a 5 x 5 matrix in 
dj and a,i + i with all other indices fixed has no symmetric traceless components. This means 
that the 14 irreps on bond < > are absent, which implies that |t/>o) is annihilated by 

the projector Hamiltonian (|3.15|) . Since the Hamiltonian (|3.15|) is positive definite, we can 
conclude that \tp ) is indeed the exact ground state. 

For a ladder with open boundary conditions, the corresponding ground state wave func- 
tion is given by 

M^a 2} ..,a N ) = T^T a \..T aN (3.25) 

This implies 4 edge states at each end of the ladder, giving rise to a 16 fold ground state 
degeneracy. 

To our knowledge, this is the first exact solution to a problem of interacting magnon and 
Cooper pairs. The ground state described by (|3.22| ) and ( [3.23D is translationally invariant 



in the ladder direction and is an 5*0(5) singlet state. It has short ranged antiferromagnetic 
and superconducting order along the ladder. Because it is a resonating state, it is hard to 
draw a simple picture for this state. The spin part of this wave function can be basically 
visualized as resonating between the states depicted in Fig. 5a and Fig. 5b. 



IV. COLLECTIVE MODES 

As discussed in Section IIIA in the -EVphase the magnon dispersion relation about q x = n 
is the same as the one hole pair dispersion relation around q x = 0. Here we examine what 
happens when the system becomes superconducting upon doping with a finite concentration 
of hole pairs. 

15 



The rung operator Li 5 (x) is equal to 



L 15 (x) = \ (N e (x) - 2) 



(4.1) 



with 



N e (x) = ]T (ct(x)c s (x) + dt(x)d s (x)) 



(4.2) 



s 



Thus Q = J2 x Li5(x) counts the number of pairs relative to half-filling and in the presense 
of a chemical potential /i, as discussed in Section II, one adds the term 



to the ladder Hamiltonian H. In the E -ph&se, as /i becomes increasingly negative, the Ei 
quintet splits with the A\Q) mode linearly decreasing its energy, the A^|Q) mode linearly 

— * 

increasing its energy and the N\Q) magnon modes remaining constant in energy, as illus- 
trated in Fig. 6a. When 2\fi\ becomes greater than the spin gap A sg , the ladder becomes 
doped with a finite density of hole pairs. As discussed in Section III, these hole pairs behave 
in strong coupling like a dilute hard core bose system with a near neighbor repulsion. In 
one dimension, hard core bosons can be treated as spinless fermions. Therefore, one can 
imagine that this band is filled with spinless fermions up to some Fermi energy, and 2|/x| 
can be physically identified with this Fermi energy. The physical origin of the increase in 
the energy to add an extra hole pair is due to the hard-core repulsion with the hole pairs in 
the condensate. 

On the other hand, the magnon band has a band minimum at momentum it and it is 
completely empty. A naive argument would suggest that one could insert a magnon simply 
by putting it at the band minimum which would cost energy A sg . However, this argument 
neglects the repulsive interaction between the magnon and the hole pairs in the condensate. 
In an 5*0(5) model, the interaction between the magnon and the hole pairs are the same 
as the mutual interaction between the hole pairs, therefore, the energy of the magnon is 
the same as the energy of a hole pair, which is 2\ji\. Thus we arrive at a simple physical 



# M = -2/jQ 



(4.3) 
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interpretation of the energy of the magnon in the superconducting state: the magnon energy 
2 1/4 1 in the superconducting state is the sum of two contributions, the rest energy A sg to 
create a magnon and the interaction energy between the magnon and the hole pair in the 
condensate. Therefore, the energy of the spin triplet momentum 7r excitation is A sg for 
\fi\ < n c and 2|/i| for > /i c , see (Fig. 6b). 

This physical interpretation was based upon a strong-coupling picture, but a more general 
result can be obtained as follows. According to Eq. ( p.24|) for an SO (5) ladder the n operators 
are exact eigen-operators 

[H + H„ 7r a ] = 2/ivr a (4.4) 

which add a pair and generate an exact excited state with momentum (jr, it) and 5 = 1. 
Now, suppose we were to start with the ground state of H with charge Q, \i/jq(Q)), and add 
a hole pair to obtain the ground state of the electron system with charge Q — 1, \ipo(Q — !))• 
The energy cost to insert the hole pair is given by the difference in the ground state energy 
of H for Q — 1 and Q electron pairs, 2|/x| = E (Q — 1) — E (Q). On the other hand, we can 
act on \ipo(Q — 1)) with the 7r|, operator and rotate the added hole pair into a magnon of 
the Q electron pair system. For the SO (5) ladder, this rotation costs no energy, so we see 
that the energy for inserting a magnon into the Q electron pair system is 2|/i|. Furthermore, 
one sees that the 7r-mode is just the natural continuation of the spin gap magnon mode |17| . 



This corresponds to the idea of the 7r-mode in the two dimensional t — J model originally 
proposed by Demler et al |6|L2"| and studied in numerical calculation by Meixner et. al. fl~0| 



and Eder et. al. 11 



The triangular relation between the 7r resonance, the magnon, and the Cooper pair can 
be illustrated by the SO (5) representation theory discussed by Eder et. al JO]. The general 
traceless symmetric tensor irreps of SO(5) are characterized by three intergers [S z ,Q,v), 
where v is related to the Casimir charge by L 2 ab = v[y + 3). This class of eigenstates of any 
SO (5) model can therefore be represented by an SO (5) pyramid, labeled by the Cartesian 
coordinates (S z , Q, v). If we take an S z = slice of the SO(5) pyramid, and plot the resulting 
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energy diagram, we obtain Fig. 7. Each box in Fig. 7 denotes the collection of all eigenstates 
with the same 5*0(5) quantum numbers. The superconducting states lie on the ridge of the 
pyramid. In the canonical ensemble, all states with the same v are strictly degenerate, but 
the lowest energy states with Av = 1 are spaced by 2\fi\ = E (Q — 1) — E (Q). Starting 
from the ground state \tp (Q)), one can create a magnon using the N a operator, or by first 
adding a hole pair using the A operator, and then acting on the resulting \ipo(Q — 1)) state 
with the 71-^ operator. (See Fig. 7). In 5*0(5) symmetric models, this triangle closes exactly, 
and the magnon mode energy is therefore predicted to be exactly 2|/x|. 

V. CONCLUSION 

We have found that a two-leg ladder with a rung interaction characterized by an onsite U 
interaction, a rung near neighbor V interaction, and a rung exchange interaction J can have 
50(5) symmetry if J = 4(U +V). Furthermore, in the E regime, the strong-coupling ground 
state is a spin-gap insulator. In this half-filled state the equal time rung magnetization and 
rung pair field correlations are identical. In addition, and of particular importance, the 
dispersion relation of a magnon rung excitation with q x measured from n is identical to the 
rung hole pair dispersion measured from q x = 0. We have also seen that when the chemical 
potential is increased such that 2\/i\ exceeds the spin gap A sg , d- wave- like hole pairs form 
a dilute hard core bose gas with a near neighbor repulsion. The spin gap magnon mode of 
the Mott insulator evolves continuously into the n resonance mode of the superconductor. 

There are also other ground states in the U — V phase diagram such as the E% regime 
which can have a CDW state or a spin gap insulating phase which when doped has s-wave 
hole pairs. In addition, the i^i-phase at half-filling corresponds to an 50 (5)-like Heisenberg 
model with ground state gaps analogous to the 5 = 1 Heisenberg model. 

A key question remains regarding the relationship of this 50(5) ladder to the more stan- 
dard Hubbard or t — J ladders. Physically if we want J and U to be positive, this requires a 
negative rung interaction V . Furthermore \ V\ must nearly balance U in order for the system 
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to be in the physical interesting regime in which J/t < 1. Therefore, the standard ladder 
models are not exactly 50 (5) symmetric in the sense defined in this paper. However, at 
half-filling, it is likely that standard ladder models flow towards a rung singlet ground state 
in the strong coupling limit. In this work we showed that such a state is not only a total 
spin singlet, but also an SO (5) singlet. Therefore, we would expect the static correlation 
to be approximately SO (5) symmetric. Recent results by Shelton and Senechal [13|, Ba- 
lents, Fisher, and Lin fl4| and Arrigoni and Hanke fl5j show that the generic interaction 
parameters of the ladder model tend to flow towards the SO (5) symmetric manifold under 
RG. However, their results were obtained in the weak coupling regime. Clearly, there remain 
various questions such as whether one will have a sufficient renormalization flow to approach 
the SO (5) regime when the physical U is of order the bandwidth. The SO (5) symmetric 
model studied in this work offer a reference point around which departures from the 50(5) 
symmetric point can be studied systematically. It would be desirable to develop a numerical 
RG analysis to study the flow around the 5*0(5) symmetric point in the strong coupling 
limit. 
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APPENDIX A: DIRAC T MATRICES AND FIERZ IDENTITY 



The general method introduced by Rabello et. al to construct 5*0(5) symmetric models 
uses the five Dirac T matrices T a (a = 1, .., 5) which satisfy the Clifford algebra, 



{r a ,r 6 } = 25 



ab 



(Al) 



Rabello et. al introduced the following explicit representation which is naturally adapted for 
discussing the unification of AF and dSC order parameters, 



(A2) 



Here a = (a x , a y , a z ) are the usual Pauli matrices and t a denotes their transposition. These 
five r a matrices form the 5 dimensional vector irreps of 5*0(5). Their commutators 





f ° 




p(2,3,4) = 


/ _ 

a 


0^ 


r 5 = 


f° 








J 















^ab 



pa yb 



(A3) 



define the 10 dimensional antisymmetric tensor irreps of 50(5). In the above representation, 
the 10 r a6 's are given explicitly by 



^15 



'-10^ 



V 

r<^ +1 ) = e ijk 

( 

V 
/ 



1 

/ 



a k 
-V fe 



p(2,3,4)l = 



— <J<Jy 



—<J y a 



{ij = 1,2,3) 



(7„ 



—ida,. 



J 
\ 



These V matrices satisfy the following commutation relations 



t <7 
-a 

*<r 



V 



pao pc 
pa& pcd 



2i(5 ac T b - 5 bc T a ) 

2t(5 ac T bd + 5 bd T ac - 5 ad T bc - 5 bc T ad ) 
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(A4) 
(A5) 



A very important property of the SO (5) Lie algebra is the pseudo- reality of its spinor 
representation. This means that there exists a matrix R with the following properties: 

R 2 = -1, i? f = iT 1 = l R=-R (A6) 
R Y a R = -*r a , R T ab R = l T ab (A7) 

The relations RT ah R~ l = — (r ab )* indicate that the spinor representation is real, and the 
antisymmetric nature of the matrix R indicates that it is pseudo-real. The R matrix plays 
a role similar to that of e a p in 5(9(3). In our representation, the R matrix takes the form 



R 



< 1^ 

v -i 0, 



(A8) 



The sixteen T A = 1, T a , T ab matrices form a complete basis in the space of 4 x 4 Hermitian 
matrices. This basis is orthonormal by virtue of the trace operation: 

Tr{T A T B ) = 45 AB (A9) 

Therefore, any 4x4 Hermitian matrix M a/ g can be expanded as 

M Q/3 = ]T X A T^ (A10) 

A 

with 

\ A = ^Tr(MT A ) (All) 

This observation can be used to derive a series of Fierz identities, relating interactions in 
the scalar, vector and tensor channels. For example, the fermion bilinear ^ a (x)^f p(y) can 
be expanded as 

= 7 D^W^d/))^ ( A12 ) 

Using this Fierz identity, one can show that the scalar, vector and the tensor interactions on 
the same rung are not independent of each other. They are instead related by the following 
equation 

(tfttf - 2) 2 = 4 - -(*t r a^2 _ I^tr^tf) 2 (A13) 
5 5 
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FIG. 1. The 16 states of a rung are laid out in their 50(5) multiplets. Here the two lines in 
the kets represent the two sites of a rung which can be in four states (empty, one electron spin up, 
spin down, or two electrons). The energies of the multiples for J = 4(U + V) are also listed. 
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FIG. 2. The half-filled strong-coupling U — V phase diagram showing the regions in which the 
single rung manifolds indicated have the lowest energy. As discussed in Section III, the Eq manifold 
is a spin-gap phase which has d x 2_ y 2-\\ke power law pair field correlations when it is doped. The E3 
manifold is divided into a CDW phase and a spin-gapped phase which exhibits s-wave like pairing 
correlations upon doping. The E\ regime is an SO (5) generalization of the Heisenberg spin one 
chain. 
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FIG. 3. Schematic illustration of (a) the Eq ground state, (b) a spin S = 1 magnon, (c) a hole 
pair. The ellipses signify a singlet state. 
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FIG. 4. Illustration of the two-degenerate CDW £"3 ground states which can occur when K3 > h. 
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FIG. 5. Schematic of two local configurations between which the system resonates in the E\ 
superspin phase. 
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(a) (b) 

FIG. 6. (a) The 50(5) vector modes versus fi. Here A + |S1) and A - add or remove a pair, 
while N\fl) represents the three 5 = 1 magnon modes. A sg = J is the spin gap. (b) The q = (ir, ir) 
magnon versus \x. These modes remain constant until \i exceeds \i c = A sg /2 and a finite pair hole 
density is formed. When this happens, the energy of the magnon- ir mode increases in energy as 

2H- 
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FIG. 7. An S z = slice of the 50(5) pyramid. Here Q = (N e - N)/2 is half of the electron 
charge measured from half-filling, and v is related to the 50(5) Casimir charge. Each box contains 
many states with the same 50(5) quantum numbers. The ground states of a given charge sector 
lie within the edge of this triangle. One can create a magnetic excitation of a charge Q state either 
by the N a operator, or by first inserting a hole pair using the A operator and then "rotate" the 
resulting state by the 7r£ operator. This triangular relationship is indicated on the figure. 
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